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INSTANTANEOUS DISTORTION IN A MACH 2.5, 40-PERCENT- 
INTERNAL-CONTRACTION INLET AND ITS EFFECT 
ON TURBOJET STALL MARGIN 
by Paul L. Burstadt and James E. Calogeras 
Lewis Research Center 

SUMMARY 

An experimental wind-tunnel investigation was conducted to determine the effects of 
time -varying total -pressure distortions produced in a supersonic inlet on a J85-GE-13 
turbojet engine. Results are presented principally in terms of "instantaneous” distortion 
amplitudes and total -pres sure contours measured through compressor stall. They indi- 
cate that although a time -averaged distortion may be far from a stall -inducing value, cor- 
responding instantaneous distortion amplitudes can approach or exceed this value. 

Seven engine stall events were studied. The time -varying data were averaged to ac- 
count for engine response characteristics so that an existing steady -state distortion pa- 
rameter and stall correlation could be applied on an instantaneous basis. When this was 
done, instantaneous distortion peaks of sufficient magnitude to cause stall were observed 
for six of the cases in the time period shortly before stall was detected. However, in- 
stantaneous distortion peaks of equal or greater magnitude were observed at earlier 
times. The fact that these earlier peaks did not cause stall indicates a lack of detailed 
understanding of the effects of distortion on compressor operation. For this reason, it 
was concluded that the instantaneous distortion approach discussed in this report might 
be used to determine a statistical probability of stall, but it does not provide an explicit 
indication of stall. 


INTRODUCTION 

Total -pressure fluctuations produced in a supersonic inlet are a major problem as- 
sociated with inlet -engine compatibility. When considered at any instant, such time- 
variant pressures can result in a large airflow distortion at the, inlet -engine interface 
relative to the time -averaged distortion. A major shar e of the fluctuations appear to 
originate from an unsteady interaction of the terminal shock with the diffuser boundary 



layer. Time-variant distortion can seriously degrade compressor stall margin. Several 
experimental investigations have already verified both the cause and effect of time var- 
iant distortion (e.g. , refs. 1 to 3). Further, such studies have qualitatively related the 
amount of stall margin degradation to the amplitude of the dynamic activity. But the pre- 
ponderance of data thus far accumulated were recorded either in an environment of arti- 
fically created distortion dynamics or in the absence of an engine. 

The present investigation was conducted in the 10- by 10 -Foot Supersonic Wind Tun- 
nel of the NASA Lewis Research Center. Its purpose was to measure time-variant dis- 
tortions produced in a supersonic inlet just before compressor stall. An axisymmetric, 
mixed -compression inlet designed for Mach 2. 5 was coupled to a J85-GE-13 turbojet en- 
gine. Steady-state and dynamic data were recorded over a range of operating conditions 
at Mach numbers of 2. 5, 2.6, and 2. 7. Results are presented principally in terms of 
’’instantaneous” distortion amplitudes and total -pressure contours measured for several 
compressor stalls. The instantaneous distortion amplitudes are compared with a steady- 
state distortion correlation to see if instantaneous distortion levels can be identified as 
the cause of stall. Compressor -face total -pressure contours are shown for each stall, 
and the relation between steady -state, instantaneous, and rms total -pressure variations 
is discussed. 


SYMBOLS 


DPCWR 

LSPR 

M 

m 

N 

N* 


NxlOO 



P 

AP 

P 

T 

W 


distortion index (defined in RESULTS AND DISCUSSION) 

loss in stall compressor pressure ratio (defined in RESULTS AND DISCUSSION) 

Mach number 

mass flow 

engine speed, rpm 

rated engine speed, 16 500 rpm 

corrected engine speed, percent 

total pressure, N/m 2 (lbf/ft 2 ) 
fluctuating component of total pressure 
static pressure, N/m 2 (lbf/ft 2 ) 
total temperature, K (°R) 
engine airflow, kg/sec (lbm/sec) 
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W corrected engine airflow, W t/<9/6, kg/sec (lbm/sec) 

cor r v 

a angle of attack, deg 

/3 spoiled sector angle, deg 

5 local corrected total pressure, P/101 325 N/m^ (P/2116 lbf/ft^) 

6 local corrected total temperature, T/288. 2 K (T/518. 7° R) 

cr standard deviation 

r data averaging period; averaging time (sec) x N/(60 sec/min), 

rotor revolutions 

Subscripts: 

0 free stream 

2 compressor -face station 

2. 1 compressor first stage exit station 

2. 2-2. 7 etc. 

3 compressor discharge station 

min, 60° minimum value in a 60° sector of station 2 
ring inner or outer annulus of station 2 
Superscript: 

_ spatial average 

APPARATUS AND PROCEDURE 
Inlet Description 

The inlet used in this investigation was designed for Mach 2. 5 and sized to match the 
airflow requirements of a J85-GE-13 engine. A photograph of the inlet terminated by a 
coldpipe -choked plug assembly is shown in figure 1 mounted in the wind tunnel. Some of 
the more pertinent details are illustrated in figure 2; more complete inlet design and 
performance details are given in reference 4. The inlet was an axisymmetric mixed- 
compression type with 40 percent of the supersonic flow area contraction provided inter- 
nally. This was accomplished with a biconic centerbody that was translated to effect in- 
let start and off -design operation. Internally, the cowl lip shock was cancelled at its 
impingement point on the centerbody. Further isentropic compression resulted in a 
throat Mach number of 1. 30. Provisions were made for both boundary -layer bleed and 
vortex generators on cowl and centerbody surfaces in the region of the geometric throat. 
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The compressor -face station was segmented by three centerbody support struts, which 
extended forward about half the length of the subsonic diffuser. Each of the three duct 
segments contained two ejector and two overboard bypasses which provided engine cool- 
ing and off -design airflow matching, respectively. Both the ejector and overboard by- 
passes were located in a common cavity, A dual -vane cascade was installed at the en- 
trance to each cavity to alleviate a resonance associated with these volumes (ref. 5). 


Engine Description 

The General Electric J85-GE-13 is an afterburning turbojet with a high thrust -to - 
weight ratio. The engine consists of an eight-stage axial -flow compressor coupled to a 
two -stage turbine. It incorporates controlled compressor interstage bleed and variable 
inlet guide vanes, a throughflow annular combustor, and an afterburner (not used in this 
test) with a variable area primary exhaust nozzle. The engine inlet diameter is 40. 9 
centimeters (16. 1 in. ). 

The exhaust nozzle area was manually controlled for this investigation. The com- 
pressor was stalled by slowly closing the nozzle while maintaining a constant engine 
speed. To avoid exceeding the turbine temperature limit during this procedure, the 
first -stage turbine nozzle area was approximately 14 percent smaller than the standard 
nozzle. This meant that at any point on the compressor map the turbine was matched to 
the compressor at a lower turbine -inlet temperature. Compressor interstage bleed, 
mechanically linked to the inlet guide vanes, was scheduled linearly from fully closed at 
90 percent of corrected engine speed to fully open at 76 percent of corrected speed. 

This was considered to be the limit of safe deviation from the normal schedule , which 
varied between 94 and 80 percent of corrected speed. This modified schedule made the 
compressor operate with less than normal stall margin. 


Instrumentation 

Steady -state and dynamic pressure instrumentation at the compressor face and com- 
pressor discharge stations is shown in figure 3(a), Only the fluctuating component of 
pressure was recorded from each of the dynamic probes. Absolute pressures were later 
obtained by adding the fluctuating component to the steady -state value provided by the 
adjacent probe. An additional correction was made at this time to account for leakage of 
the high -pass filters, which had caused a small zero -shift of the fluctuating signals. 
Frequency response of the dynamic probes at the compressor face was flat td about 
2000 hertz. Response of the compressor discharge dynamic probes and the interstage 
static probes shown in figure 3(b) was flat to about 300 hertz. 
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Acquisition and Reduction of Data 


A block diagram of the analog data acquisition system is shown in figure 4(a). The 
fluctuating signal from each pressure transducer was recorded on FM -multiplexed tape 
at 152.4 centimeters per second (60 in. /sec). This speed provided recording capability 
of up to 4000 hertz. 

Figure 4(b) is a block diagram of the data reduction system for the 30 compressor - 
face total pressures. The FM tape was played back at 9.525 centimeters per second 
(3. 75 in. /sec), or 1/16 of the original recording speed. Signals from the tape were de- 
multiplexed and then passed through a bank of linear -phase amplifier -filters . Gain set - 
tings were selected to provide the correct ranges for the digitizer, and bandwidth set- 
tings were selected to provide low -pass filtering at an effective frequency of 1600 hertz 
(i.e. , 16 times the bandwidth setting). The data in a 2 25 -millisecond period (about 215 
msec before stall) were then digitized at an effective rate of 8000 points per second per 
channel and put on tape. This digitizing rate provided five points per cycle of the highest 
frequency signal that passed through the amplifier -filters. 

Calculations were made for each time slice generated by the digitizer. Reduced data 
were then provided in the form of microfilm plots and/or printed output. 

The rms amplitude of the fluctuating component of total pressure was determined for 
each of the compressor -face dynamic probes. This was done on-line by passing the dy- 
namic signals through low -pass filters with 1000-hertz corner frequencies, and then 
measuring them with true rms (trms) meters. The output of each meter was then re- 
corded on the steady -state data system after passing through another low -pass filter with 
a corner frequency of 1/24 hertz. 

Dynamic data from the compressor interstage and discharge measuring stations 
were not processed through the reduction system shown in figure 4(b) but were reduced 
to the form of oscillograph traces, using the linear -phase amplifier -filters. 


Data Bank 

A description of the set of stall events recorded in this program is shown in figure 5. 
For many of the stalls the event recorded on FM tape and the corresponding steady -state 
data scan were recorded at slightly different operating conditions. The test procedure 
was to establish an engine operating condition in the vicinity of the compressor -stall 
point. The FM tape recorder was started and a steady -state data scan was made. Then, 
while the FM tape recorder was still running, the compressor pressure ratio was in- 
creased slightly and the throttle adjusted to maintain a constant rotor speed (see Engine 
Description). If stall did not occur, this procedure was repeated until the stall event 
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was recorded on the FM tape. In some cases the stall occurred after the engine had sta- 
bilized at an equilibrium operating condition. These were called ’’drift” stalls. 

Seventeen of the 29 stalls recorded were of the ’’drift” type, and in 11 cases a 
steady-state data scan had been completed immediately before stall. These 11 cases 
provided the most accurate data for determining distortion values, and for this reason 
they were the only cases examined in detail. Two of these 11 cases had relatively low 
levels of both steady -state and dynamic, (^ p rms /^2’ distortion. As a result, they 
showed very little loss in stall margin. Two other cases showed significant losses in 
stall margin, but the dynamic level was low, and the loss was attributed to their respec- 
tive steady -state distortion levels. So of the 29 stall events recorded on tape, only 
seven were determined to be suitable for instantaneous distortion analysis. 


RESULTS AND DISCUSSION 
Steady -State Compressor Performance 

Compressor performance with both undistorted and hub-radial distorted inflows is 
presented in figure 6. The stall compressor -pressure ratios from these figures are 
shown in figure 7 as a function of corrected engine speed. These data were recorded in 
an altitude test facility as part of a program to determine the effect of steady -state total- 
pressure distortions on a J85-GE-13 turbojet (see ref, 6). The same engine was used in 
the wind tunnel test program discussed herein, so that correlations from the steady-state 
program could be applied without considering engine -to -engine variations. 

Figure 7 is the basis for determining the loss in stall compressor pressure ratio 
(LSPR) suffered by the engine at any stall point. The compressor map for hub -radial 
distortion is substantially different from either the undistorted map or those resulting 
from circumferentially distorted inflows shown in reference 6. For this reason it was 
decided to base LSPR on the dashed curve of figure 7, whenever steady -state hub-radial 
distortion was present, and on the solid curve for all other cases. Then LSPR is defined 
as: 


LSPR = 


P,/P. 


1 - 


distorted stall 


P 3^2 


ref 



= constant 



where (P 3 /P 2 ) ref is from figure 7. 
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Steady-State Compressor Sensitivity to Stall 


The sensitivity of the compressor to various types of steady -state distortion patterns 
is presented in figure 8. These results are from the screen-induced distortion program 
reported in reference 6. Loss in stall compressor pressure ratio is plotted as a function 
of a distortion index based on a critical angle of spoiled flow. This index, DPCWR, was 
developed from the test results of reference 9 and is determined from the two inner rings 
and two outer rings of compressor -face total -pressure probes shown in figure 3(a). 

A simple radial averaging procedure is used to convert the 12 pressures in the two 
inner rings into an equivalent set of six pressures, each representing a 60° sector. The 
same procedure is then applied to the two outer rings. The distortion is then calculated 
for each (inner and outer) equivalent ring as 


distortion = 



> min, 60° j 

P / . 

ring 



where 

P ring average pressure 

P min 60° pressure in ring 

^ring 60° times number of contiguous pressures in ring less than P 

The maximum of these two distortion values is called DPCWR. 

The data in figure 8 represent stall points for each of the patterns shown. The dis- 
tortion index takes no direct account of pure radial distortions. But the effect of hub- 
radial distortions is considered in the definition of LSPR, as mentioned previously. 
Midspan radial distortions had little effect on compressor performance, and they were 
neglected in the definition of DPCWR. Tip-radial distortions did affect compressor per- 
formance, but the effect was not considered since an accurate correlation was not 
available. 

In figure 8 the rms error, or one standard deviation a, of the DPCWR values with 
respect to the faired curve was determined to be 0. 01. For random data scatter the 
±0. 01 band would include 68 percent of the data points. So when a value of DPCWR is 
determined from a given value of LSPR, an error band of at least ±0. 01 must also be 
considered. This error band can be quite large in terms of percent for small values of 
DPCWR. However, we will be considering only the absolute value of the error band 
when discussing the accuracy of a correlation between instantaneous distortion levels and 
compressor stall. 
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Description of the Seven Drift Stalls 


Inlet -engine operating conditions are listed in table I for the seven drift stalls to be 
analyzed. The predicted critical distortion in table I (also referred to as the critical 
level) is the value of DPCWR from figure 8 corresponding to the value of LSPR, which 
accounts for hub -radial distortion as discussed previously. It is the value of distortion 
index predicted to be necessary to cause the measured LSPR and has an associated error 
band of ±0. 01, as discussed earlier. 

The average rms amplitude of the fluctuating component of total pressure is also 
shown in table I. Each of the values recorded from the trms meters (see Acquisition and 
Reduction of Data) was divided by Pg, and then an average of all 30 ratios was computed. 

Steady -state operating points for the seven drift stalls are shown plotted on a com- 
pressor map in figure 9. The speed characteristics for undistorted inflow are shown as 
dashed lines. Corrected speed values for the stall points are listed in table I. 

A series of compressor -face total -pressure contours for each of the seven drift 
stalls is shown in figure 10. The contours represent the steady -state , dynamic (i.e. , 
rms), and instantaneous total pressures. They were made by fitting a surface through 
the 30 individual pressure values and then assigning shadings to different ranges. Al- 
though the aft section of the subsonic diffuser was segmented by three centerbody support 
struts, they were ignored when the contours were made. The instantaneous total - 
pressure contours correspond to high levels of instantaneous distortion and will be dis- 
cussed in conjunction with subsequent figures. They are included in figure 10 for ease of 
comparison with the steady-state contours. 

The descriptions of the steady -state distortion patterns given in table I can be seen 
from examination of figure 10. Contours of the fluctuating component of total pressure 
(rms) indicate that the highest level of dynamic activity (dark shading) generally occurred 
in a band between the high and low total -pres.sure regions. This band corresponds to a 
region with large spatial total -pressure gradients. These regions of high dynamic activ- 
ity are strongly related to the instantaneous distortion amplitudes that preceeded com- 
pressor stall. 


Determination of Time of Stall 

To evaluate the results of the instantaneous distortion analysis, it was necessary to 
have an accurate indication of the time when stall occurred inside the compressor. The 
signals from the compressor static -pres sure transducers (see fig. 3(b)) were put out on 
strip charts and examined. For all seven drift stalls the row of instrumentation at the 
290° location showed the earliest indication of stall and these are the charts shown in 
figure 11. Since the origin of a stall 'one in a compressor stage causes an abrupt loss 
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in the airflow pumping of that stage, a compression wave is sent forward and an expan- 
sion wave is sent rearward. Previous investigations with the J85 (ref. 7) have indi- 
cated that stall originates in a circumferential sector of the compressor. The stalled 
sector then increases in circumferential extent as it rotates in the same direction as the 
rotor, but at about half rotor speed. Results indicate that, in general, the stall zone 
makes one or two complete revolutions before terminal flow breakdown occurs, causing 
the rapid rise in inlet total-pressure known as the hammershock. This means that the 
hammer shock cannot be used as an accurate indication of the stall time, since the period 
between the first indication of stall and the hammershock is variable. 

Because there were only two axial rows of instrumentation inside the compressor, 
it was not possible to determine exactly when stall originated inside the compressor. 
Although the first indication was seen at the 290° location, a stall zone could originate 
just past the 45° location but not be detected until it rotated past the instrumentation at 
290°. Since the stall zone rotated at about half rotor speed (8 msec/rev) there could be 
as much as a 6 -millisecond delay before it was detected at the 290° location. This pos- 
sible time delay will be important when discussing the instantaneous analysis in later 
figures. 

Figure 11(a) shows that for reading 141, stall was first observed inside the com- 
pressor at 43 milliseconds. (The zero of the time scales corresponds to the start of in- 
stantaneous distortion calculations. ) Static pressure rose at the exit of stages 1 and 2, 
and pressure dropped at the exit of stages 4 through 7. (The initial disturbance is not 
clearly shown at the stage 3 exit. ) This indicates a stall zone which originated in stage 3 
or 4. Complete flow breakdown occurred at about 50 milliseconds. 

Charts for the other six readings are shown in the remaining sections of figure 11. 

In some cases the initial indication of stall is not as obvious as the disturbance caused 
by the second passage of the stalled zone. Figure 11(e) is an example of this situation. 
The first and third stages show small disturbances at about 36 milliseconds. Since 
stronger disturbances were observed about 7 to 8 milliseconds later (the time required 
for one rotation of the stalled zone) , it was concluded that the initial small disturbances 
were a valid indication of stall. 


Effect of Data Averaging on Time -Varying Distortion 

In figure 12 the steady -state distortion index (DPCWR) is applied on a time -variant 
basis for a typical wind tunnel stall point (reading 148). Although 225 -milliseconds of 
digitized data were available before each stall, only the last 40 -milliseconds (10 rotor 
revolutions) before stall were analyzed to reduce the amount of data processing and to 
allow more detailed examination of the time period near stall. 
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The 50 -millisecond time span used in figure 12 represents a period just before and 
including the complete breakdown of flow in the compressor (terminal stall) and is iden- 
tical to the time span used in figure 11(b) for reading 148. To compose this figure, the 
distortion computations were made and values plotted every 0. 125 milliseconds. The 
critical level and its error band, along with the stall time and steady -state distortion 
level, are shown on the plots. Figure 12 shows increasing values of the time increment 
over which compressor -face pressures were averaged before computing distortion 
values. When time -averaging distortion data, the averaging procedure should be applied 
directly to the pressure data rather than to the time -varying distortion index. It is clear 
from figure 12(a) that if the time varying index is averaged, it will not necessarily ap- 
proach the steady -state value as a limit. 

The averaging time of 0.07 rotor revolution (approximately 0. 25 msec) shown in 
figure 12(a) represents the effective averaging time of the predigitizer analog filters 
shown in figure 4. (More details on analog averaging are given in ref. 8. ) The averag- 
ing times used in figures 12(b) through (e) represent the effects of a digital averaging 
procedure as well as the effects of the analog filters. For these cases the digital (i. e. , 
arithmetic) average pressure value for each probe was determined from the digitized 
values encompassed by the averager "window” (of width r) and was assigned the value of 
time corresponding to the center of the ’’window. ” The averager window was then moved 
forward by an increment of 0. 125 millisecond, and the process was repeated. 

As ejected, increasing the averaging time r reduces the level of the instantaneous 
distortion peaks, smooths the variation of distortion with time, and reduces the mean 
value of the fluctuations. The mean value of DPCWR approaches the steady -state level 
as r approaches infinity. A more subtle effect of increasing r is that the relative 
magnitudes of the distortion parameter peaks are not maintained. For example, the 
highest distortion peak before stall in figure 12(a) occurs at 18 milliseconds, but in fig- 
ure 12(c) the highest peak occurs at 5 milliseconds. This effect is important when the 
response of the compressor to variations in distortion is considered. 

Since the instantaneous distortion approach is simply the application of a steady - 
state distortion parameter on a quasi -steady -state basis, the most ideal representation 
of this approach would be to associate compressor stall with a single distortion peak, 
which has the following characteristics: 

(1) It is the highest distortion peak that occurs in the span of time beginning with 
equilibrium engine operation and ending with the first indication of stall inside the 
compressor. 

(2) The time increment between the highest distortion peak and the first indication 
of the stall should be no longer than the transport time required for the distortion to 
reach the stalling stage, plus the time between the origin of stall and the first indication 
by the interstage instrumentation, plus the time delay in blade -element response (i.e. , 
flow separation). 
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(3) The amplitude of the highest instantaneous distortion peak should equal or exceed 
the critical level of distortion, as determined by the steady -state correlation of figure 8. 

Reference 9 mentions a delay of less than 0. 5 millisecond for blade element re- 
sponse. For the results presented in this report a period of 1 millisecond was used to 
represent the combination of transport time and blade-element response. As mentioned 
previously, the time increment between the origin of stall and the first indication of that 
stall could be as long as 6 milliseconds. In this investigation the first indication of stall 
was usually observed while the stalled zone was still rather small. It was not until some 
later time, when the stalled zone had grown in size, that the resulting compression wave 
(hammershock) affected the compressor -face total pressures. And, since the hammer - 
shock did not affect all pressure rakes at the same time (since the stalled zone was still 
rotating), a large increase in distortion amplitude usually resulted. This corresponds to 
a time of about 47 milliseconds in figure 12. 

To meet the third criterion, it is necessary to time average the fluctuating compo- 
nent of compressor -face total pressure for a period that represents the dynamic re- 
sponse of the engine. The compressor would then be expected to respond to the resulting 
instantaneous distortions in a quasi -steady -state manner. The ratio of maximum instan- 
taneous to critical distortion can be used to evaluate the effect of averaging time. When 
the ratio is 1.0, the highest level of distortion equals the critical level, and it is as- 
sumed that the steady-state stall criterion (fig. 8) can be applied. 

Figure 13 is a plot of the ratio of maximum instantaneous to critical distortion as a 
function of averaging time for each of the seven drift stalls. The basic aspects of 
DPCWR have been programmed on an analog computer (reL 10), and the results show 
that high levels of the parameter exist at many times before stall. But it was felt 
that the number of high distortion levels that existed in the 40 millisecond -period before 
stall were sufficient to allow a valid determination of the averaging time required to 
bring the peak values down near the critical level. 

The square symbols (reading 148) in figure 13 represent the data from figure 12 and 
the steady-state distortion. Plots similar to figure 12 were available for all readings 
but are not shown in the report. For all seven readings the ratio of steady -state to crit- 
ical distortion is represented by a solid symbol in figure 13, plotted at the value of 
r = 00 . The set of solid symbols shows clearly that the steady -state levels were not suf- 
ficient to cause stall, based on the correlation of figure 8. 

Examination of figure 13 shows that five of the seven readings had their desired av- 
eraging times in the range of 0. 2 to 0. 75 rotor revolution. This amount of averaging is 
assumed to be an indication of the engine’s dynamic response. It appears that reading 
141 would have reached the desired level of peak distortion for a value of r equal to 
about 1. 5 rotor revolutions. This larger averaging time appeared to result from the 
high level of dynamic activity ((APj^g/P)^ = 0. 085) present at the compressor face. On 
the other hand, reading 164 had no instantaneous distortion values above the critical 



level during the time period examined. The best single value of r for all seven read- 
ings appears to be 0.5 rotor revolution. 

Evaluation of Time -Varying Distortion for Seven Drift Stalls 

Figure 14 shows one instantaneous -distortion plot for each of the drift stalls. The 
plots were selected based on the desired values of r determined from the curves in fig- 
ure 13. Because the data were only processed for five values of t, the plots are for r 
close to the desired value. 

Time varying distortion for reading 141 is shown in figure 14(a). The critical level 
and its associated error band are shown in the figure. The band encompasses a large 
segment of the data partly because the value of r = 1. 18 for this plot is somewhat less 
than the desired value of about 1. 5 rotor revolutions. But only two of the distortion peaks 
shown can be related to the stall event. 

A large value of distortion occurred just as stall was observed inside the compres- 
sor, at 43 milliseconds. But this distortion value is partially dependent on events that 
occurred after the stall was observed. In this case the averager ’’window” was about 4 
milliseconds (1. 18 rotor revolutions) wide. The distortion computed from averaged 
pressure values is plotted at the time corresponding to the midpoint of the window, so the 
peak at 43 milliseconds was influenced by pressures recorded between 41 and 45 milli- 
seconds. As mentioned before, the effects of stall (hammershock over -pressures) did 
not usually appear at the compressor face until a few milliseconds after stall was ob- 
served inside the compressor. But for reading 141 the hammershock appeared at the 
compressor face at the same time stall was first observed inside the compressor. This 
could be due to the possible time delay, mentioned earlier, between the origin of stall 
and the initial observation of stall inside the compressor. So, instead of only being in- 
fluenced by pressures recorded after stall was observed, the distortion peak at 43 milli- 
seconds was also influenced by the hammershock. It’s clear that results must be ex- 
amined carefully when averaging periods exceed the time between the distortion peak and 
the observation of stall inside the compressor . 

After considering the 1 -millisecond period for transport time and blade response and 
the possible 6 -millisecond delay between the origin and observation of stall, it is as- 
sumed that a distortion peak should be seen in figure 14(a) between about 36 and 42 
milliseconds. 

There is a distortion peak equal to the critical level at 36 milliseconds and another 
at about 39 milliseconds. Although these peaks do not represent the maximum value of 
distortion calculated before stall, they are of sufficient magnitude to cause stall (con- 
sidering the ±0. 01 band) and they occur at times consistent with the origin of stall inside 
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the compressor. The large number of peaks near the critical level may mean that this 
index gives a statistical, rather than an explicit, indication of stall. 

Figure 10(a) shows instantaneous compressor -face total -pressure contours that cor- 
respond to two of the distortion peaks in figure 14(a). These peaks occurred too early to 
have caused stall. However , these contours are representative of patterns observed for 
many other distortion peaks. They show that the peak levels of instantaneous distortion 
are caused by pressure contours which are more circumferential in nature than the 
steady-state pattern. They look like hub -radial patterns that have been stretched along a 
diameter. It seems clear that this stretching is related to the band of high dynamic ac- 
tivity in the rms contour (also fig. 10(a)). 

The distortion time -history for reading 148 is shown in figure 14(b) for an averaging 
time of 0. 31 rotor revolution, which is almost enough to bring the maximum distortion 
down to the critical level. Two distortion peaks (at about 34 and 37 msec) of amplitude 
sufficient to cause stall occurred during the 7 -millisecond period before stall was ob- 
served inside the compressor. A somewhat smaller distortion peak occurred at 33. 5 
milliseconds, and it also could have been the cause of stall. But it occurred very early 
in the allowable time period and if a slightly longer averaging time had been used (to 
reach the desired case of maximum == critical distortion) the peak would probably fall out- 
side of the ±0. 01 band about the critical level. (Compare figs. 12(c) and (d). ) So it 
seems reasonable to assume that the peak at 34 or 37 milliseconds was responsible for 
the stall. It might seem obvious to assume that the later peak (37 msec) caused the stall. 
But it is possible that the peak occurred after stall had already originated inside the com- 
pressor. Stall did not affect the compressor face (via the hammer shock) until much 
later. 

The instantaneous total -pressure contours corresponding to the distortion peaks near 
34 and 37 milliseconds are shown in figure 10(b), When these are compared with the 
steady -state contour also shown in figure 10(b), it is clear that the instantaneous patterns 
have much larger and less symmetric regions of low pressure. These low-pressure 
bulges also occurred at other circumferential locations at other times. Examination of 
the rms total -pressure contour in figure 10(b) again indicates that the large pressure 
fluctuations represented by the dark band are the probable source of the differences be- 
tween steady-state and instantaneous patterns. 

Reading 154 also had a steady -state hub -radial distortion pattern, and its instanta- 
neous distortion plot is shown in figure 14(c). The distortion peak at about 32 millisec- 
onds seems to be the obvious cause of stall. Plots for larger values of t (not shown) 
indicate that this peak would remain the point of interest for a slightly larger t neces- 
sary to bring it down to the critical level. 

The maximum value of distortion for reading 154 occurred about 0. 5 millisecond be- 
fore the stall was observed inside the compressor. -Since this would be the minimum 



allowable time for transport and blade response, it is more likely that stall was caused 
by an earlier value of distortion that existed on the ’’uphill” side of the peak. 

The instantaneous compressor face contour for the maximum distortion point at 
about 32 milliseconds is given in figure 10(c) along with tne steady-state and rms pres- 
sure contours. Although the steady -state and. rms contours are very similar to those of 
reading 148, the instantaneous contours are slightly different. But the instantaneous 
contours in figure 10(c) still appears to be basically a hub -radial pattern that has been 
stretched by the action of the band of high dynamics. The high distortion value was 
caused by the combination of high and low pressures in the outer annulus . 

A distortion time -history and the various compressor -face pressure contours for 
reading 162 are shown in figures 14(d) and 10(d), respectively. For an averaging time 
of 0. 59 rotor revolution, three distortion peaks had magnitudes close to the critical 
level. But based on the time stall was observed, it appears that stall was caused by a 
distortion that existed on the uphill side of the peak at 38 milliseconds. 

Since the inlet was at a = 5° for reading 162, the steady -state distortion pattern 
was primarily circumferential and this was also the case for the instantaneous patterns. 
The highest level of dynamic activity (i. e. , rms level) occurred in the high pressure 
gradient region and resulted in an expansion and contraction of the high-pressure region. 
The instantaneous contour in figure 10(d) shows an expanded high-pressure region (with 
respect to the steady -state pattern) along with some regions below the minimum steady- 
state levels. A 16 -millimeter film was made using a complete series of instantaneous 
contours from reading 162, and it showed that the major variation between contours was 
due to the changing extent of the high-pressure region along with occasional appearances 
of regions less than the minimum steady -state pressure levels. 

The expansion and contraction of the high-pressure region in the film indicated that 
when the inlet is at angle -of -attack, the line between high- and low-pressure regions 
sweeps back and forth across the compressor face probes in a rather cyclic manner. A 
similar film was made for reading 148 (a = 0°), and the variation of the instantaneous 
pressure patterns was observed to be of a more random nature. 

The operating conditions for reading 164 were very similar to those of 162, as can 
be seen from table I. The corrected speed was lower for reading 164, so the inlet t er- 
minal shock was farther forward. This gave a higher Pg/Pg and lower (AP rms /P) 2 - 
But the LSPR, and therefore the critical distortion, was about the same for both 
readings. 

The instantaneous distortion plot for reading 164 is shown in figure 14(e), and did 
not involve any digital averaging. But even for this lowest averaging time, the distortion 
never reached the region within ±0.01 of the critical level. 

Figure 10(e) shows two instantaneous pressure contours for reading 164, and their 
relation to the steady-state and rms contours (same fig. ) is similar to that observed for 
reading 162 in figure 10(d). 



The ±0.01 error band about the critical level only represents a ±1 ct error band for 
the steady-state distortion correlation of figure 8. But even if a ±3 a (±0. 03) band was 
placed in figure 14(e), it would not include any distortion peaks. And since a ±3cr band 
would be expected to account for essentially all random data scatter in figure 8, it ap- 
pears that the present distortion analysis is not adequate to explain this stall event. 

For reading 103, vortex generators were installed on the inlet centerbody, but the 
steady-state distortion pattern was still basically hub-radial in nature. Figure 14(f) 
shows the distortion time -history for this reading using an averaging period of about one- 
half rotor revolution. When the previously discussed time allowances are considered, it 
seems probable that the distortion pattern at about 31 milliseconds was the cause of stall. 

Examination of figure 10(f) shows that the instantaneous pressure contour is substan- 
tially different from the steady-state pattern and shows much greater circumferential 
distortion. The rms contour of figure 10(f) again shows the peak fluctuations in the re- 
gion of high-pressure gradient, but the contour is less symmetric than those observed 
for the other cases of steady -state hub -radial distortion. It is possible that this could be 
a result of the vortex generators. The low-pressure region at the outer edge of the in- 
stantaneous contour corresponds to a region of low dynamic activity. But it is also near 
a low-pressure bulge in the steady-state pattern, which may explain why this low instan- 
taneous pressure was observed in a region of low dynamic activity. 

Figure 14(g) is the type of plot that clearly displays the basic aspect of the instan- 
taneous distortion concept; a single peak in the distortion index which greatly exceeds all 
others and occurs just before engine stall. This figure is the instantaneous distortion 
variation for reading 261 (a = 6°), and the very sharp definition of the peak near 34 milli- 
seconds is caused by an approximation in the computer program used to calculate the dis- 
tortion index. 

As mentioned in the definition of DPCWR, the extent functions /3 r j n g only varies in 
increments of 60°. This means that significant errors can result for certain orientations 
of regions which are below the average pressure. For reading 261 the value of 
varied between 120° and 180° as the less -than -average pressure region moved back and 
forth across the probes. (Nearly all the distortion values plotted in fig. 14(g) correspond 
to the outer two rings of probes. ) This stepwise fluctuation of the extent function is the 
cause of the sharply defined distortion peaks near 30 and 32 milliseconds , the large peak 
near 34 milliseconds, and some other step changes that occur at earlier times. 

If the computer program had been designed to calculate the extent function in a more 
accurate and continuously varying manner, the results would have been quite similar to 
figure 14(g). However, the variation of distortion index with time would have been much 
smoother and the absolute values would be somewhat different. But' the general features 
of figure 14(g) are considered to be valid. 

The various total -pressure contours for reading 261 are shown in figure 10(g). The 
instantaneous contour corresponds to the highest distortion peak and shows an extent of 
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pressure below average of between 120° and 180°. Again, the high rms pressure levels 
occur in the region of highest steady -state total-pressure gradient. 

For all of the drift stalls except reading 164, an instantaneous distortion peak of 
magnitude sufficient to cause stall was observed less than 7 milliseconds before the 
observed stall time. The 7 -millisecond period was the maximum time increment that 
could be expected to exist between occurrence of a stall inducing pattern and the first in- 
dication of stall. The magnitude of distortion sufficient to cause stall was based on the 
steady-state distortion correlation (fig. 8) and its associated data scatter. To account 
for the data scatter, all values of distortion determined from figure 8 were assumed to 
include a ±la error band of ±0. 01. More recent tests of other J85-13 engines in an alti- 
tude chamber (ref. 11) indicate that a ±lcr error band much larger than ±0. 01 would be 
necessary to cover the scatter of data obtained from several J85-GE-13 engines. 


SUMMARY OF RESULTS 

An experimental investigation was conducted to determine the effect of time -varying 
total -pressure distortions on a J85-GE-13 engine. The engine was run behind a mixed- 
compression inlet in the Lewis 10- by 10 -Foot Supersonic Wind Tunnel. Data were re- 
corded from high-response total pressure probes at the compressor face. A previously 
developed steady -state distortion index was then applied to the digitized pressure data on 
an instantaneous (quasi -steady) basis, and an attempt was made to relate high levels of 
instantaneous distortion to compressor stall. The following conclusions were drawn: 

1. Although the steady -state distortion level may be far from the predicted stall - 
inducing value, corresponding instantaneous distortion amplitudes can approach or ex- 
ceed this value. 

2. When the inlet is at zero angle -of -attack, severe total -pres sure fluctuations can 
result in an instantaneous distortion pattern that is quite different from the steady -state 
pattern. For angle -of -attack cases the instantaneous distortion patterns are similar to 
the steady-state circumferential pattern. 

3. The region of highest total -pressure fluctuations was in a band located between 
the high and low (steady -state) pressure recovery zones. A great deal of the variation 
between instantaneous pressure contours seemed to be the result of these regions of high 
total -pressure fluctuations. 

4. It is necessary to have an accurate indication of the time stall first occurs inside 
the compressor to help identify the distortion that caused stall. 

5. The appearance of the hammershock at the compressor face is not an accurate 
indication of the time stall occurs because the time delay between the first indication of 
stall inside the compressor and the appearance of the hammershock is not the same for 
all engine stalls. 
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6. For this J85 engine the time -varying distortion data need to be averaged for a 
time period on the order of 0. 5 rotor revolution before the steady-state stall criterion 
can be applied on an instantaneous basis. 

7. After applying the proper averaging period and when the data scatter of the 
steady-state correlation was considered, an instantaneous distortion peak of amplitude 
sufficient to cause stall was observed for six of the seven drift stalls less than 7 milli- 
seconds before stall was first observed inside the compressor. 

8. The existence of numerous ’’false indications” of stall indicate that the present 
parameter and instantaneous distortion approach might provide a type of statistical prob- 
ability, instead of a positive indication, that stall will occur. This is caused by a lack 
of detailed understanding of the effects of steady-state (and, therefore, instantaneous) 
distortion on compressor operation, and the current necessity of using a primarily em- 
pirical approach to correlate distortion effects. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, February 5, 1974, 

501-24. 
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pressure 

□ Static 

o Total 

Open symbols - steady state 
Solid symbols - dynamic 



Compressor face Compressor discharge 

(looking downstream) 


(a) Compressor face and compressor discharge measuring stations. 


no 



(b) Compressor interstage measuring stations. 


Figure 3. - Steady-state and dynamic pressure instrumentation. 




Digitized data 


Total number 
of stalls re- 
corded on FM- 
multiplex tape 


Number of 
'drift' stalls 


Number of 
'drift' stalls 
recorded on 
steady-state 
system 


Number of 
digitized stalls 
with appreci- 
able loss in 
stall margin 


Number of 
digitized stalls 
with significant 
levels of dynam- 
ic distortion 



Figure 5. - Data bank resulting from 60/40 inlet-J85 engine compatibility study. 
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Compressor pressure ratio, P3/P2 




Figure 6. - Compressor performance. 



Figure 7 . - Variation of stall compressor pressure ratio with engine speed. 
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Loss in stall compressor pressure ratio, LSPR 



0 © © 





Compressor pressure ratio, P 2 /Pg 

0.900 TO 0.850 0.800 TO 0.750 ** 0.700 TO 0.660 


0.850 TO 0.800 m 0.750 TO 0.700 * 0.650 'AND B&OH 

(a-1) Steady-state total -pressure contour. Aver- 
age compressor pressure ratio, P 2 /Pg, 0.761 



Compressor pressure ratio, P^Pg 

0.900 TO 0.850 " ! 0. 800 TO 0.750 * 0.700 TO 0.850 

0.850 TO 0.800 ^O.TSOTO 0.700 "0.650 AND BELOW 


(a-3) Instantaneous total-pressure contour. 
Time, 20.50 milliseconds; averaging time, 
1. 18 rotor revolutions. 


rms pressure ratio, (A P rrns /P ^2 

88 0.1 to TO 0.095 *-0.080 TO 0.065 

*0.125 TO 0.110 “ 0.095 TO 0.080 0.065 TC 0.050 


(a-2) rms t otal -press ure contour. Dynamic dis- 
tortion, {AP rms /P) 2 , 0.085. 



Compressor pressure ratio, P 2 /P 0 

0.900 TO 0.850 K! 0.800 TO 0.750 aB 0.700 TO 0.650 

^ 0.850 TO 0.800 88 0.750 TO 0.700 *0.650 At© BELON 


(a -4) instantaneous total-pressure contour. 
Time, 33.75 milliseconds; averaging time, 
1. 18 rotor revolutions. 


(a) Reading 141. 

Figure 10. - Steady-state, dynamic and instantaneous total -pressure contours. 



Compressor pressure ratio, P 2 / Pq 

0.900 TO 0.050 0.800 TO 0.750 * 0.700 TO 0.650 , 

0.850 TO 0.800 88 0.750 TO 0.700 


(b-1) Steady-state total -pressure contour. Average 
compressor pressure ratio, P 2 /Pq, 0.788. 



Compressor pressure ratio, P 2 /P 0 

0.900 TO 0.850 555 0.806 TO 0.750 * 0,700 TO 0.650 

■* 0.850 TO 0.800 88 0.750 TO 0.700 

(b-3) Instantaneous total -pressure contour. 
Time, 33.875 milliseconds; averaging time, 

0. 31 rotor revolution. 



rms pressure ratio, (AP rms /P) 2 


88 0.080 TO 0,070 ^0.060 TO 0.050 

*0,090 TO 0.080 *0.070 TO .0. 060 0.050 TO 0.0*0 

(b-2) rms t otal -press ure contour. Dynamics dis- 
tortion, {AP rms /P) 2 , 0.066. 



Compressor pressure ratio, P^Pg 

0.900 TO 0.850 “ 0.800 TO 0.750 *0.700 TO 0,650 

m 0.850 TO 0,800 88 0.750 TO 0.700 

{b-4) Instantaneous total -pressure contour. 
Time, 37.375 milliseconds; averaging time, 
0.31 rotor revolution. 


(b) Reading 148. 
Figure 10. - Continued. 
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Compressor pressure ratio, P^Pg rms pressure ratio, (AP rms /P )2 


6.900 TO 0.850 “ 6.800 TO 0.750 * 0.700 TO 0.650 *0.090 Mb ABOVE *0.080 TO 0.070 ^ 0.060 TO 0.050 

^naRftTO 0.800 *0.750 TO 0.700 * 0.090 TO 0.080 ss 0.070 TO 0.060 0.050 TO 0.040 

(c-1) Steady-state total-pressure contour. Aver- (c-2) rms t otal-press ure contour. Dynamic dis- 
age compressor pressure ratio, P 2 / Pq. 0.799. tortion, (AP rms /P) 2 , 0.070. 



Compressor pressure ratio, P 2 / Pq 


0.900 TO 0.050 * 0.800 TO 0.750 *0.700 TO 0.650 

® 0.050 TO 0.800 *0.750 TO 0.700 

(c-3) Instantaneous total-pressure contour. 
Time, 32.375 milliseconds; averaging time, 
0. 33 rotor revolution. 


(c) Reading 154. 
Figure 10. - Continued. 
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Compressor pressure ratio, P 2 /Pq 

0.9*0 TO 0.880 "*0.820 TO 0.760 *0.700 TO 0.6*0 

23 0.880 TO 0.820 *®0.760 TO 0.700 

(d-1) Steady-state total -pressure contour. Aver- 
age compressor pressure ratio, P 2 /Pg, 0.767. 


rms pressure ratio, (AP rms /P) 2 

^O. 060 TO 0.050 0.0*0 TO 0.030 

m 0.070 TO 0.060 0.050 TO 0.0*0 0.030 TO 0.020 

(d-2) rms total-pressure contour. Dynamic dis- 
tortion, (AP rms /P) 2 , 0.046. 



Compressor pressure ratio, P 2 /Pq 

C.900 TO 0.850 =* 0.800 TO 0.750 *0.700 TO 0.650 

0.850 TO 0.800 **0.750 TO 0.700 *0.650 At® BELO* 

(d-3) Instantaneous total pressure contour. 
Time, 37.875 milliseconds; averaging time, 
0.59 rotor revolution. 

(d) Reading 162. 

Figure 10. - Continued. 
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Compressor pressure ratio, P 2 /Pg 

0.940 TO 0.080 383 0.820 TO 0.760 

^ Q.880 TO 0.820 m 0.760 TO 0.700 


(e-1) Steady-state total pressure contour. Aver- 
age compressor pressure ratio, P 2 /Pg, 0*844. 



Compressor pressure ratio, P 2 /Pq 


rms pressure ratio, (AP rms /P) 2 

*0.030 AND ABOVE ^0.02510 0.020 " 0.015 TC 0.310 

* 0.030 TO 0.025 5K 0.020 TO 0.015 :;i .0.fl10 ?C 0.305 


(e-2) rms t otal -press ure contour. Dynamic dis- 
tortion, (AP rms /P) 2 , 0.019. 



Compressor pressure ratio, P 2 /Pg 


0.900 TO 0.850 555 0.800 TO 0.750 

P-s 0.050 TO 0.800 88 0.750 TO 0.700 

(e-3) Instantaneous total -pressure contour. 
Time, 35.25 milliseconds; averaging time, 
0.066 rotor revolution. 


■ : 0.900 TO 0.850 0.800 TO 0.750 

0.850 TO 0.800 88 0.750 TO 0.700 

(e-4) Instantaneous total -pressure contour. 
Time, 36.125 milliseconds; averaging time, 
0.066 rotor revolution. 


(e) Reading 164. 


Figure 10. - Continued. 
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Compressor pressure ratio, P 2 /Pg 


rms pressure ratio, (AP rms /P) 2 


“ 0.050 TO 0.0*0 0.070 TO 0.0S0 

0-0*a TC 0.330 “ 0.060 TC 0.050 

(f-2) rms to tal -press ure contour. Dynamic dis- 
tortion, (AP rms /P) 2 , 0.050. 



0.780 TO 0.760 “ ! 0,7*0TO 0.720 **0.700 TO 0,680 

531 0.760 TO 0.740 *0.720 TO 0.700 *0.680 AND BELOH 

(f-1) Steady -state total -pressure contour. Aver- 
age compressor pressure ratio, P 2 / Pq- °- 736 - 


Compressor pressure ratio, P^Pq 


* 0.820 TO 0.760 *0.700 70 0.6*0 

*0.760 70 0.700 

(f-3) Instantaneous total -pressure contour. 
Time, 32.125 milliseconds; averaging time, 
0.55 rotor revolution. 

(f) Reading 103, 

Figure 10. - Continued. 


30 


4 



i i iitffl 



Compressor pressure ratio, P 2 /P 0 rms pressure ratio, (A P rms /P )2 


0,900 TO 0.850 0.800 TO 0.750 * 0.700 TO 0,650 88 0.015 TO 0.055 “ 0,025 TO 0.015 

m 0.850 TO 0.800 *5 0.750 TO 0.700 *0.650 AND BEL OH “O.OKTO 0.025 0.015 TC 0.005 

(g-1) Steady-state total -pressure contour. Aver- (g-2) rms t otal-press ure contour. Dynamic dis- 
age compressor pressure ratio, P/Pq, 0.783. tortion, (AP rrns /P> 2 , 0.021, 



Compressor pressure ratio, P 2 /Pq 

■•0.900 TO 0.850 28 0.800 TO 0.750 * 0.700 TO 0.650 

m 0.850 TO 0.800 * 0.750 TO 0.700 * 0.650 W© BQ.OH 


(g-3) Instantaneous total-pressure contour. 
Time, 33.250 milliseconds; averaging time, 
0. 19 rotor revolution. 


(g) Reading 26. 
Figure 10. - Concluded. 
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Static pressures 
at exit of com- 



0 10 20 30 40 50 


(d) Reading 162 


Static pressures 
at exit of com- 
pressor stages 



0 10 20 30 40 50 


Static pressures 
at exits of com- 
pressor stages 



Time, msec 


(e) Reading 164. 


(f) Reading 103. 


Figure 11. - Continued. 


34 




Static pressures 
at exits of com- 
pressor stages 



0 10 20 30 40 50 


Time, msec 
.(g) Reading 261. 
Figure 11. - Concluded. 
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Distortion index, DPCWR 




0 10 20 30 40 50 



(b) Reading 148; averaging time, 0. 19 rotor revolution. 



0 10 20 30 40 50 


(c) Reading 148; averaging time, 0.31 rotor revolution, (d) Reading 148; averaging time, 0.55 rotor revolution. 



Time, msec 


le) Reading 148; averaging time, 1.03 rotor revolutions. 
Figure 12. - Effect of averaging on time-varying distortion. 







Ratio of maximum instantaneous to critical distortion 




Distortion index, DPCWR 



(a) Reading 141; angle of attack, 0 ; corrected engine (b) Reading 148; angle of attack, 0 ; corrected engine 
speed, 98.82 percent; compressor pressure ratio at speed, 86.88 percent; compressor pressure ratio, 

stall, 7.045; averaging time, 1.18 rotor revolution. 4.601; averaging time, 0.31 rotor revolution. 



0 10 20 30 40 50 0 10 20 30 40 50 


Time, msec 

(c) Reading 154; angle of attack, 0 ; corrected engine (d) Reading 162; angle of attack, 5°; corrected engine 
speed, 92.25 percent; compressor pressure ratio at speed, 92.89 percent; compressor pressure ratio at 

stall, 5.240; averaging time, 0.33 rotor revolution. stall, 5.288; averaging time, 0.59 rotor revolution. 

Figure 14. - Time-varying distortion for wind tunnel stall points. 
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0 10 20 30 40 30 0 10 20 30 40 50 

Time, msec 


(e) Reading 164; angle of attack, 5°; corrected engine 
speed, 86. 89 percent; compressor pressure ratio at 
stall, 4.388; averaging time, 0.07 rotor revolution. 


(f) Reading 103; angle of attack, 0; corrected engine 
speed, 86. 42 percent; compressor pressure ratio at 
stall, 4.684; averaging time, 0.55 rotor revolution. 



(g) Reading 261; angle of attack, 6°; corrected engine 
speed, 87. 00 percent; compressor pressure ratio, 
4.355; averaging time, 0.19 rotor revolution. 


Figure 14. - Concluded. 


MASA-Lanelev, 1974 
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